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Abstract 

We analyse time series of CDS spreads for a set of major US and European institutions on a pe- 
riod overlapping the recent financial crisis. We extend the existing methodology of e-drawdowns to 
the one of joint e-drawups, in order to estimate the conditional probabilities of abrupt co-movements 
among spreads. We correct for randomness and for finite size effects and we find significant prob- 
ability of joint drawups for certain pairs of CDS. We also find significant probability of trend rein- 
forcement, i.e. drawups in a given CDS followed by drawups in the same CDS. Finally, we take the 
matrix of probability of joint drawups as an estimate of the network of financial dependencies among 
institutions. We then carry out a network analysis that provides insights into the role of systemically 
important financial institutions. 



Introduction 



Within the field of complex networks (Caldarelli 2007 1, the investigation of financial networks is cur- 



rently one of the emerging avenues (Schweitzer etal.| 2009 ), also in view of the on-going global financial 



crisis. Financial contagion on networks (Batt iston et al.|[2012b| ) differs in some important respects from 
the well-known processes of epidemics spreading ( [Barrat et al.| 2008 Pastor-Satorras and Vespignani , 



2001 ). It also bears similarities to epidemics spreading, among which, the fact that the topological struc- 
ture of the network plays a crucial role in the collective dynamics and therefore in emergence of systemic 
risk. There is a body of work on networks reconstructed from correlations among equity prices or return 



time series (Bonanno et al. 2003 Garasetal. 2008 Kullmann et al. 2002). The analysis of the minimum 
spanning tree provides insights into the classification of stocks and the level of correlation depending on 
the market phase. Correlation analysis suffers, however, from some important limitations, the main one 
being that zero correlation between two series does not imply that they are independent (only the inverse 
is true, see more details in section 4 in the SI). To overcome these limitations, here we utilise a method 
based on the detection of joint £-drawup, which allows us to estimate the probability that two series 
exhibit a co-movement. Moreover, in contrast to equity, CDS prices reflect the default probability of the 
reference entity and thus the network constructed from CDS prices are more relevant in studying the 
propagation of default risk. 

Our approach can be applied to construct networks of dependencies in other financial markets. In general, 
it applies to all domains of networks in which links are, for any reason, unobservable but the dynamics of 
the nodes reflect the dependency structure. To summarise, the contributions of the paper are the following. 
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First, we build on the e-drawup method ( Sornett e and Zhou] 2006 1 to estimate the probability of joint 
e-drawups, which are essentially a particular type of co-movements across time series. Based on this, 
we estimate the level of the so-called interdependence and trend reinforcement in the system across 
different phases of the market. In addition, we find them to be relevant indicators for the emergence 
of systemic risk, as was also predicted by previous theoretical work. Second, we construct a network 
of interdependencies among institutions and we introduce two novel centrality measures that allow the 
identification of systemically important nodes in the network. Our approach enables the disentanglement 
of a structure that is, only apparently, very homogenous. It also allows us to track the role of nodes 
evolving in time. 



1 Results 

Credit Default Swaps (CDS's) are financial derivatives instruments in which the seller provides the buyer 
protection against a credit event of a reference entity (see section 2 of the SI). We have analysed the 
time series data of CDS prices, or spreads, of top US and European financial institutions from January 
2002 'till December 2011. We observe three different phases in the CDS prices, corresponding to the 
periods before, during, and after the credit crisis of 2008: a normal phase (period 1); a volatile and 
upwards trending phase (period 2), and a volatile and downwards trending phase (period 3) (see figure 
[T]). Accordingly, we divide the data into three parts. For each institution's time series we detect the so 



called e-drawup's. An £-drawup is an extension of the notion of an e-drawdown (Johansen (2003 1). It 
refers to a persistent upward movement in a time series until a peak has been reached, after which the 
time series declines (or, has a "correction") by more than an amplitude e (see figure[5jb)). Since the CDS 
spread represents the cost of insurance, an e-drawup signifies an increase in the default probability of 
that institution, as perceived by the market. When market participants buy and sell insurance on each 
other, their financial performances can become interdependent (see section 3 of the SI). 

Whilst interdependence can be seen as a form of risk diversification which decreases individual risk, 



previous work (Battist on et al.|[2012a i has demonstrated that high interdependence leads, instead, to 



higher systemic risk when coupled to a so-called trend reinforcement. 
Interdependence and trend reinforcement. 

In our context, trend reinforcement refers to the tendency of an e-drawup to be followed by another e- 
drawup in the same time series. Interdependence refers, in contrast, to the existence of co-movements 
between two different time series (i.e. an e-drawup followed by another one in a different time series). 
Here, we take the frequency of e-drawup 's in i as an estimate of the probability Pj that security / has 
a e-drawup. Similarly, for the frequency of joint e-drawup's we estimate Py, i.e. probability that j ex- 
periences a e-drawup given that i experiences a e-drawup. We also account for a time lag t = 0, 1, 2, 3 
days between the drawup's. The expected probability of joint drawup's in the case of two statistically 
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Figure 1 : Time series of credit default swaps throughout the credit crisis. A plot of the CDS spread 
time series covering the financial crisis of 2008. The data ranges from January 2002 to December 201 1. 
We can observe three market phases. Most CDS spreads peak around March 2009. The CDS prices are 
quoted in basis points (bp). 

independent time series is P ; j = PiPj. Therefore, we take as an estimate of interdependence between two 
financial institutions, the deviation from such a case, i.e. Wij = Ptj - PiPj- Finally, in order to account 
for finite size effects, we consider only the values of Wij that cannot be rejected based on a permutation 
test and we reset to zero all the other values (see Methods). We show the distribution of Wj/s and W;,-'s 
(Fig.[T^), b). Notice that, because of the time lag, Wjj is not a symmetric matrix. A positive value of Wij 
does not prove any causality between the movements of i and /, but measures the dependence of j due 
from i, or equivalently, the impact of i on j. Accordingly, the histograms (Fig.[T^), b) only show positive 
values of W,/s and Wa's. 

Analogously, we take as an estimate of the trend reinforcement for institution i the deviation: Wu = 
Pu - PjPj and we treat it as above. 
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We find a significant level of trend reinforcement (i.e. which cannot be accounted by chance or finite 
size) in about 50%, 72% and 80% of the nodes (respectively in period 1, 2 and 3). This is reflected in 
the histograms of Wu (Fig. [jj>). As we can see, the curve for period 2 and 3 is mostly above the one for 
period 1. This means that the number of nodes with a significant level of trend reinforcement increases 
when the market moves from the first phase to more volatile phases. We also find a significant level of 
interdependence in 54%, 78% and 77% of pairs of nodes in period 1, 2, 3, respectively. The histograms 
of Wij (see Fig. []J>) show that periods 2 and 3 are characterised by higher frequencies. In fact, 20% of 
pairs in period 2 and 19% pairs of nodes in period 3 exhibit values of Wij greater than the mean plus one 
standard deviation of period 1. Moreover, while in period 1, nearly all values of Wij are smaller than 0.5, 
in period 2 and 3 there is a tail extending up to 1 . 

These findings show that interdependence and trend reinforcement are indeed present in an important 
market such as the one for CDS's. Moreover, trend reinforcement increases from period 1 to period 2, 
and even more so does interdependence. Before drawing the implications of these findings in terms of 
systemic risk, we proceed to a network analysis of the structure of interdependencies among institutions. 
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Figure 2: a)The distribution of non-zero values of interdependence Wij across the three periods. The 
counts in periods 2 and 3 are higher than in period 1 . In addition, during periods 2 and 3 the distribu- 
tions of Wij have longer tails compared to period 1 . b) The distribution of non-zero values of trend 
reinforcement Wu across the three periods. 



Network analysis and systemic risk. The CDS market can be naturally mapped into a directed and 
weighted network in which nodes represent institutions and edges represent interdependencies among 
institutions. More precisely, whenever Wij > 0, we assign a weighted edge with value Wij from institution 
i to j. Since Wij is the probability that conditional to i experiencing a draw-up, j also experiences a draw- 
up with a time lag t, it follows that the stronger the edge, the stronger the impact that i has on j. 
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When looking at the properties of connectedness of the network, we find a significant number of dis- 
connected nodes in all three periods (81, 39, 39, respectively). Remarkably, the rest of the nodes form 
only one strongly-connected component (LSCC, see Methods) encompassing, respectively, 95, 137, 137 
nodes in period 1, 2 and 3. The density of links (i.e. the number of links over the number of possible 
links) in the LSCC is high in all the three periods: 0.98,0.97, 0.97. This is reflected also in the average 
out degree in the LSCC's across the three periods, which is 90 ± 8, 129 + 11, and 123 ± 19. Finally, the 
average path length within the LSCC's is 1.04, 1.05 and 1.2, meaning that almost all the nodes in the 
LSCC are first neighbours to each other. In such a structure, each node has a direct impact on all the other 
nodes, and each of these has a further impact on all the others. Intuitively, this finding suggests that the 
financial distress at one node in the SCC can quickly propagate to all the other nodes in the LSCC and 
keeps reverberating through the many connections. 

Indeed, recent works on systemic risk in financial networks have shown that the number of links play an 
ambiguous role. Few links are functional to diversify the individual risk. However, too many links gen- 
erate systemic risk. This holds in presence of mechanisms that either amplify the distress (such as in the 
case of contagion), or simply increase the persistency of the distress in time, such as trend reinforcement 
(Battiston et al.||2012a[ ). As we have seen, the CDS market exhibits a core of more than 100 nodes, that 



is almost a fully connected graph (i.e. with maximal degree) and where in many cases links represent 
strong interdependencies. According to the theoretical results mentioned earlier, in such a situation, even 
small levels of amplification can make the whole system very unstable. Notice that moving from period 
1 to period 2, the values of most CDS's raised dramatically, in many cases by one order of magnitude 
(see Fig.[T]). Looking at individual institutions it is clear that their risk of default had become very high. 

However, previous findings suggest that, in such a situation, the default of a few institutions would have 
triggered a systemic default. Indeed, it is generally thought that without the massive intervention of the 
Federal Reserve (FED) through various emergency programs that lasted from the fall of 2008 until the 



summer of 2009, (Bloomberg-News 201 1 



Therefore, monitoring trend reinforcement and interdependence together with the individual riskiness of 
the participants seems to provide a valuable assessment of the level of systemic risk in a market. 

Centrality. In order to gain insights into the systemic importance of specific nodes in the network, we 
proceed to investigate their centrality. We obviously focus only on the strongly connected component 
since the other nodes are isolated. 

The out- and in-degree of a node are the simplest measures of centrality that hold a valuable interpretation 
here: A high out-degree represents the ability of a node to affect many neighbours when it experiences 
a draw-up; a high in-degree corresponds to a node being affected by many nodes. Since the network 
is almost a complete graph, based on the out-degree, all nodes are equally systemically important and 
equally affected by the others. Interestingly, even the out- and in-degree strength of a node (i.e. the sum 
of the weights in the outgoing and incoming links), yield similar conclusions. 



5, 



15 



Eidgenbssische Technische Hochschuie Zurich Rahul Kaushik and Stefano Battiston: 

Swiss Federal Institute of Technology Zurich Kaushik and Battiston, CDS Networks 

http : //www . sg . ethz . ch 



As an alternative approach, based on the notion of feedback centrality, for each node i we introduce a 
novel measure of systemic importance, called impacting centrality and denoted as b,. The measure takes 
into account, in a recursive way, the fact that a node is more systemically important if it impacts many 
systemically important nodes (see Methods). Symmetrically, we also introduce a measure of systemic 
vulnerability, called impacted centrality of a node i and denoted as c,-. The latter measure captures the 
idea that a node is more vulnerable if it has strong dependencies from many nodes which are in turn 
heavily impacted by the whole network - Recall that we are only considering upward movements in 
CDS spreads, i.e. an increase in the fragility of the underlying reference firm. 

In both cases, the values are normalised between and 1. Remarkably, in stark contrast to in- and out- 
degree, the values of these two centrality measures are broadly spread across the range [0, 1] (see Fig.[TJ 
also SI, Fig. 6a & 7a). 

If we focus on the ratio r, = ^, in the scatter plot of Fig. jlj it is possible to identify three regions (above, 
between and below the dotted lines), corresponding to three different roles. Nodes are located in the top 
region if they have a value r, > 3/2, meaning that they impact the network 1.5 times more than they are 
impacted by it. Symmetrically, nodes are in the bottom region if r, < 2/3. Finally, nodes that appear in 
the middle region are those that impact and are impacted by the network in a comparable manner. 

A scatter plot of impacting vs impacted centrality in each period is shown in Fig. [TJ where the size of 
each node reflects the average (relative) debt held by a financial institution during that period. According 
to the above classification, in period 1 most institutions are located in the middle region, implying that 
the most vulnerable nodes are also the most systemically important. Moreover, the nodes with the largest 
debt also happen to be the most systemically important and vulnerable, as they are located in the top 
right corner. 

In period 2 and 3 most institutions progressively move to the top and the bottom regions. For instance, 
in period 3 (red bubbles) we can observe one cluster nodes with large debt in the top centre corner of 
the plot. This reflects that these nodes have a high impacting centrality, but only an intermediate level of 
impacted centrality. A possible explanation is that these institutions managed to reduce their exposure to 
their counterparties from period 2 to period 3. There is also a cluster of nodes with large debts located in 
the bottom centre of the plot. These nodes possibly correspond to firms that not only managed to reduce 
their exposure, but also became less systemically important. 

We observe that there are many nodes in the network that not only have a high impacting centrality, but 
also a high impacted centrality. In addition, by adjusting the size of each node based on the average debt 
held by a financial institution a clearer picture emerges. 

Unlike period 2, financial institutions with the largest debt's were in the top right of Fig. [T] The top right 
corner of Fig.[T]represents firms that are just as vulnerable as they impact the others in the network. Since, 
the size of the debt of such institutions is the largest, a small percentage change in their debt would cause 
a large change in the distress of the debt issued by other firms in the network. We find that in contrast 
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to the evolution of CDS spreads in period 1, the impacting and impacted centralities along- with average 
debt of the institutions implied ripe conditions for greater levels of systemic risk. We find that firms 
with high average debt that are mostly being impacted are important from a systemic risk perspective. 
However from a financial regulator's perspective it is also important to track firms, based on their size, 
impacting, and impacted centrality, that have the ability to cause distress to a large subset of the market. 

From a systemic risk perspective it is essential to study nodes that are prone to distress; however, from a 
policymakers perspective it is also vital to monitor nodes that are not only prone to distress, but that also 
have a high impacting centrality as distress in such nodes could lead to a systemic collapse. 

Bow-tie extraction. In order to emphasise the 3 different roles suggested by this finding, we carry out 
the following filtering of the links. In each period, for nodes located in the top region, we remove all their 
incoming links. Symmetrically, for those in the bottom region, we remove all the outgoing links. Since 
the initial network is strongly connected, in this way. We, by construction, obtain a bow-tie structure (see 
Methods). The position of a node in the bow-tie is related to its systemic importance. Indeed, the IN, 
SCC and OUT component of the bow-tie correspond to the top, middle and bottom regions of Fig. [T] 
respectively (e.g. the nodes in the IN are those that impact the network more than they are impacted). 

Note that the bow-tie structure is constructed based on the choice of impacting-impacted centrality, i.e. 
nodes with r, > 3/2 are in the IN, nodes with 2/3 < r ; - < 3/2 are in the SCC, and nodes with r,- > 2/3 
are in the OUT. In fact, for any 8 > 0, where 8 € (0, 1). The lines 1-8 and 1+8 would separate the 
nodes into three regions. Thus, the choice of 8 is based on the level of impacting-impacted centrality that 
is of interest, see SI for more visualisations. In addition, if W is a directed SCC, and one truncates all 
incoming links of nodes with r, < 1 - 8, and all outgoing links for nodes with r,- > 1 + 8. Then, it is not 
always the case that the filtered W has a non-trivial SCC, see SI for more details. 

We then introduce a novel method for the visualisation of the bow-tie (Fig.Q. 

This enables the representation, at the same time, of a network structure, the position of the nodes in the 
various component of the bow-tie, as well as their level of impacting centrality. In the diagram [4j the 
circle represents the SCC, the top (bottom) section correspond to the IN (OUT). E.g. within the SCC, 
more central nodes are located towards the centre of the circle. The colour code and the size of the dots 
also reflect their centrality, such that the red and large dots are the most central (see caption of Fig. 
[4]). This visualisation allows to track how individual institutions become more or less central, or if they 
changed role across periods (see SI, Fig. 6 & 7 b). 

In period 1, most of the nodes of the bow-tie are in the SCC (85), with 4 and 6 in the IN and OUT 
respectively. Moreover, most nodes in the centre of the SCC are banks and investment banks, while 
insurance and real estate companies tend to be in the periphery of the SCC (see SI, Fig. 6 b). 

This implies that in the normal phase most of the nodes impact the network, and are impacted by the 
network in a comparable manner. In period 2, the bow-tie grows overall, but the SCC (97 nodes) grows 
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Figure 3: Scatter plot of impacting versus impacted centrality. Each institution in the CDS market 
is represented by three dots depending on the period (blue, green, red refers to period 1, 2, 3, respec- 
tively). The size of each node is determined by the leverage of a financial institution. It can be seen that, 
while in period 1 most institutions are located between the two dotted lines, in period 2 and 3 many 
of them move to the top and bottom region. This means that ratio between the two centrality measures 
varies with the market phase. Few institutions of interest are labelled. For example, Bank of America 
(BOFA) remains in the same region across the three periods. With reference to the subsequent bow-tie 
construction used in Fig. [4] The scatter plot is divided into 3 regions. Nodes in the region above the line 
n > 3/2 correspond to the IN. Nodes in the region 2/3 < n < 3/2 correspond to the SCC. Nodes in the 
region r,- < 2/3 correspond to the OUT. 



proportionally less than IN (19 nodes) and OUT (22 nodes). Of the 81 nodes that were disconnected in 
period 1, twenty seven migrated to the SCC (Fig. [4]). 
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Figure 4: The network of the CDS reference entities from period 2. Each of the nodes represents a 
financial institution. Outgoing links from nodes that are in the top, or the IN of the bow-tie structure 
represent the estimated potential impact of a financial institution to its neighbours (see Methods). The 
nodes in the SCC are placed within a circle of radius one and centred at the origin. The distance of 
each node from the centre is 1 -Impacting centrality. The angle increases linearly from to 2n. Thus, 
the closer a node is to the centre the higher is impacted-impacting centrality. Similarly, nodes in the 
OUT and IN are placed between angles n/2- 5tt/8 and 3n/2 - 13/r/8 respectively. In addition, nodes in 
the OUT and IN are placed with an offset of 1 . 1 from the origin. With the bow-tie representation we 
are able to visually compare the centrality of a node i with node/ Also, with this visualisation we are 
able to extract a network of nodes that mostly impact the others, nodes that impact just as much as they 
get impacted, and nodes that only get impacted by other nodes in the network. The size and the colour 
of the node reflects impacted-impacting centrality of a node (nodes with larger impacted-impacting 
centrality are in red). The colour assigned to links is based on where the links point to in the network. 
Links originating from IN to the SCC are in bright blue. Links originating in the SCC to nodes in the 
SCC are in green. Links that are originating in the SCC to the OUT are dull blue grey colour. 

In period 3, the size of the bow-tie remains unchanged, but the SCC shrinks by about a 50% (from 97 
to 47 nodes), as a result of a migration to the IN (37 nodes) and mostly to the OUT (53 nodes) (see SI, 
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Fig. 7 b). In particular, the nodes with high impacting centrality are now all located in the IN and not, 
anymore, in the SCC. 

One should not forget that the original network is a strongly connected graph and the bow-tie is obtained 
with a filtering. Therefore, it is not the case that the nodes in the IN are not connected among each 
other. This means that in case a few nodes would have defaulted, the others would still have been heavily 
affected. However, the observed migration of nodes implies that, compared to the normal period, there 
has been an increasing polarisation between nodes (IN) that predominantly impact the network, and 
nodes (OUT) that predominantly are impacted by the network. 

The above analysis of impacting and impacted centralities, and the bow-tie extraction allows us to move 
from an initial picture in which all nodes seemed to be equally important from the point of view of 
systemic risk, to a much more refined picture. In terms of systemic impact, we can now focus on a small 
subset of the nodes, viz. nodes that have a high impacting centrality and are located in the centre of 
the SCC, or in the top part of the IN. This finding is corroborated by anecdotal evidence -(see, e.g., 



FCIC 2011 )- about the role of important actors of the credit crisis of 2008 (see SI). Conversely, the 



impacting centrality allows also to identify nodes that suffer the most from an impact originating from 
the others. Remarkably, there is no evidence of one or two nodes dominating the others in terms of 
systemic importance. In contrast, we see that in each period a set of about top 19 nodes have similar 
values of centrality. 



2 Discussion 

We have analysed the e drawup's in the CDS's time series for the top US and EU institutions throughout 
the last 10 years. By measuring the frequency of joint drawup's in pairs of CDS time series we have 
estimated the level of interdependence and trend reinforcement in the market. According to previous 
theoretical works on financial networks, the interplay of these two mechanisms is deeply linked to the 
emergence of systemic risk. We have found very significant levels of both, and an increase of both in acute 
phases of the crisis. The result suggests that high interdependence and trend reinforcement together with 
high level of individual riskiness are a good indicators of the level of systemic risk. Indeed, when CDS 
spreads where at their peak in 2008, implying high risk of individual default, movements in the spread 
of a few institutions were very likely to be followed by movements in another and also in the same 
institutions. This means that the default of a few players was likely to trigger a financial melt-down. 

Furthermore, we have carried out what to our knowledge is the first study of the complex network of 
CDS interdependencies. In order to investigate the systemic importance of individual nodes, we have 
introduced two novel measures. The impacting centrality captures, in a recursive way, how much a node 
impacts the network. Symmetrically, the impacted centrality captures how much a node is impacted by 
the network. These two measures enable the extraction of a bow-tie structure from the initial network 
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and to clarify the role of the nodes. We have found that in all phases of the market the top institutions by 
systemic importance are not just one or two, but about 19. 

They all have comparable level of systemic importance and they all are very tightly interconnected. 

The specific findings of this analysis are relevant to the broad audience interested in the issue of systemic 
risk and systemically important financial institutions, including policy makers. Moreover, our approach 
is very general and applies to any set of time series associated to units that operate in interaction. In 
particular, it is of interest for those cases where the direct interaction between units is not observable and 
the dependence has to be inferred from the dynamics. In this respect, our paper contributes to a stream of 



work on the observability and the reconstruction of complex networks (Eldawlatly et al. 2010),(Baca- 
Lopez et al.|[2009l ,( |Pfitzmann and Waidner[ [T987) . 



3 Methods 



Data description The data consists of CDS spreads of single name entities denominated in US dollars 
and the Euro encompassing a total of 176 firms. 




Figure 5: a) Illustration of the e-drawup methodology. The represent local extrema that were not 
detected as drawups. The red-dots represent local maxima that were picked up as candidates for a e- 
drawup. The green-dots represent the local minima that were picked up for a e-drawup. Compare the 
difference between the maxima and minima on days 5 and 6 respectively with e. Since, e > than the 
difference, we iterate to the next set of local maxima and minima on days 7 and 8, keeping day 3 as the 
day from when we count a e-drawup b) The plot highlights the e-drawup methodology applied to the 
time series of American International Group (AIG) and Merrill Lynch (MER) 
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We have focused our attention on three different periods corresponding to different phases, see [T] of the 
market: (1) April 2003 - May 2006 (representative of a normal phase); (2) May 2006 - March 2009 
(volatile with an upwards trend); (3) March 2009 - December 2011 (volatile with a downwards trend 
market scenario). The motivation to do a period-wise analysis is to extract clear snapshots of a network 
structure before, during and after the crisis of 2008. This data window covers a period of 3624 trading 
days. 

Data consistency The data was acquired via a subscription to Bloomberg. The CDS data in general is 
provided to Bloomberg by various quote providers. Bloomberg has negotiated contracts with various 
quote providers to display data on the Bloomberg terminal; however, not all quote providers relax the 
constraint of Bloomberg being allowed to distribute this data to the end user. Thus, we see certain breaks 
in our time series, as on those days either the CDS was not traded (see SI), or is quoted by a quote 
provider that does not allow its data to be downloaded via Bloomberg. To circumvent this issue, we carry 
forward the last traded price till a time that a new price has been reported. 

e-drawups We compute the e-drawup's in each of the time series using the following algorithm, see 
figure [5] Suppose, we start our analysis of a e-drawup from the first green point from the left in figure 
[5J3). We follow the following steps to compute £-drawup's: (1) We compute the local variation in the 
time series for the last ten days, call it e. (2) compute local extrema. (3) Goto first local minima, call it 
s-drawup can didate (on day three, figure[5^). (4) Iterate to the set of local maxima and minima (occurring 
on days 5 and 6 respectively, see figure[5j)). Compute the difference between the maxima and the minima, 
call it correction amplitude {correction amplitude refers to the decline in price followed after an increase 
in price). (5) Update e. (6) If correction amplitude > e. Then, we record the e-drawup on the day it occurs 
(day 7 in figure [5Jb)). And, we update s-drawup ca ndidate- Otherwise, we iterate to the next minima and 
goto the succeeding maxima and repeat steps above. 

In this way, we are able to compute the e-drawup 's in the the time series data. Once we have computed 
the e-drawup 's in the time series we map it into the vector Vi of node i. Where the length of the vector is 
the same as the length of the time series, T. Also, Vj(f) — 1 if there was a drawup on day t, otherwise we 
put it to zero. After proceeding with this methodology, we extract N vectors of e-drawup's, where N is 
the number nodes in the analysis into a matrix of drawup's, w,, where i e {1, ...,N). 

Co-movements Once we have computed the matrix, w,, we proceed to computing a matrix of co- 
movements. We implement the following algorithm to compute co-movements: (1) We select an arbi- 
trary but fixed node i. (2) We loop from t = 1 till t — T and compare each Wi(t) with all wj(t) where 
je {1 , N}. Compare with all w j(t + t) where ; € {1, N] and r € {0, 1,2,3}. (3) Ifw,-(0 - 1 and 
Wj(t) = 1 , then count = count + 1 . If If w ,■(/) - 1 and wj(t + t) = 1 , then counf - counf + 1 , and then end 
the iteration. (4) Update matrix of joint drawup's, i.e. D, 7 = DT. = C0U j f . Repeat steps above till all 
the i, j pairs have been accounted for. We then filter each of the e-drawup's based on a permutation test. 
In order to estimate the probability that a node j experiences a e-drawup in the next three days given that 
node i experiences a e-drawup today, we sum the filtered matrices, i.e. Xr-n ^ r ■ Notice that the quantity 

1 yJ I J 
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represented by Wjj is not a measure of causality, as each Wij represents, in a probabilistic framework, the 
probability of co-movements and not that node i causes a e-drawup in node j. 

The matrix of co-movements now carry probabilities of joint drawup's in the time series for all pairs. 
The matrix of co-movements, constructed as above, is square; however, D, y , DJ. are not symmetric. 

Randomness. To compute each element of the matrix Wij we subtract (element wise) from matrices 
Dij and D„ * Djj and D r H * D T .. respectively. Alternatively, by the following procedure as above, we 
essentially conduct a test of independence, i.e. Wij = Pij - P%Pj- In addition, by construction Wjj is 
adjusted for random occurrences. 

Finite size. In order to clean the matrix of drawup's, Wu. We permute the matrix indices of where the 
drawup's occur in w, . We permute each pair of vectors w\ and wj hundred times. After each permutation 
we re-compute Wf° ntrol (i, j) hundred times, following the procedure (see Methods), and extract a single 

Wj/"""" 1 from hundred realisations at the 95% confidence interval. We then filter each element in Wu 

~ control 

using our results from Wij . We then proceed to interpreting this filtered weighted matrix as an 
adjacency list of a directed network. 

Impacted and Impacted centrality In line with the notion of feedback centrality (e.g., PageRank see 
SI) we define the following measures. 

bi^WijVjW+fi^Wijbj (1) 

j j 

j j 

where Wij = YiW-' ~ ^' anc ^ P = ^-85 is a dampening factor. The matrix Wij is row stochastic 



and thus equations [2] and[l] 
b = (I- W) _1 Wfc. Equation 



can be analytically solved to yield the solution, c = (I - W') : W'v and 
T]can be interpreted as the importance of a firm i based on how much it gets 



impacted by the other firms in the system. Equation [2]can be interpreted as the centrality of a firm based 
on how it affects other members in the network. Where Wy is the transpose of the adjacency matrix that 
we derive from the e drawup procedure and after having cleaned for randomness and finite size. 

The impacting centrality is analogous to: cumulatively how much of the distress in the network is due to 
node i. Or, the loss experienced by the other nodes by the failure of node i. Suppose there are ,/V shocks 
that are experienced by node i that are propagated to the other nodes in the network based on Wjj. Thus, 
by counting the number of times a shock from i passes through the other nodes in the network represents 
the impact of node i on the network. Similarly, by counting the number of shocks in the entire system 
that propagate through node i represent the impacted centrality of node i. 

The impacting centrality can be viewed as the cumulative sum of the shocks induced on all the nodes in 
the network cumulatively over time due to node i. OR, the expected no. of shocks that are propagated 
into the system given that I have a shock impacting centrality. 
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impacted centrality is the cumulative fraction of the shocks induced in me by a set of shocks on all the 
other nodes in the network over time. OR, the expected no. of shocks that are propagated through me by 
all the other nodes in the system. 

Bow-tie structure A bow-tie network is essentially a directed network consisting of four main bodies, 
namely a Strongly Connected Component (SCC); OUT: set of all nodes that can be reached, directly or 
indirectly, from the SCC; IN: the set of all nodes that reach the SCC directly or indirectly. The fourth 
and last component of the bow-tie structure, Tubes and Tendrils (TT) represent the set of all nodes that 
are not a part of the SCC; however, a node in the TT can either be reached from the IN and/or OUT. If 
r,- > 3/2 then we remove all the incoming links. This means that all nodes that exhibit r, > 3/2 only have 
outgoing links. Similarly, for nodes that exhibit r,- < 2/3 have all the outgoing links removed. And, the 
remaining nodes, i.e. 2/3 < ri < 3/2 are in the SCC and retain both the incoming and outgoing links. 
Clearly, this does not imply that nodes in the IN, or OUT do not have links among themselves; but, we 
simply disregard them by construction. This merely highlights the nodes ability to impact its neighbours 
more than it gets affected and vice versa. Thus, all nodes in the IN point to nodes in the SCS, and nodes 
in OUT get pointed to by the nodes in the SCC. 
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